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ALBORAN  SEA  MODELING 


1.  INTRODUCTION 

The  circulation  of  the  Alboran  Sea  is  dependent  on  both 
Atlantic  and  Mediterranean  waters.  Atlantic  water  flows  through  the 
narrow  (20  km  wide)  and  shallow  (300  m  deep)  Strait  of  Gibraltar  into 
the  Alboran  Sea  forming  a  150-200  m  deep  surface  layer  (Ovchinnikov, 
1966;  Lanoix,  1974;  Katz,  1972).  Mediterranean  water  enters  the  Alboran 
at  its  open  eastern  boundary  and  flows  slowly  westward  in  the  form  of  an 
intermediate  and  lower  layer.  It  has  been  suggested  that  even  the 
deepest  water  can  exit  through  the  Strait  (Stommel  et  al.,  1973; 

Gascard,  1982).  The  large  volume  transport  of  inflowing  Atlantic  water, 

1  to  2  X  106  m3/sec,  (Lacombe,  1971;  Bethoux,  1979;  Lacombe,  1982) 

retains  its  identity  as  a  narrow  (30  km  wide)  jet  with  initial  speeds  of 
100  cm/sec  near  the  Strait  (Lacombe,  1961;  Peluchon  and  Donguy,  1962; 
Grousson  and  Faroux,  1963;  Lanoix,  1974;  Cheney,  1977;  Petit  et  al., 
1978;  and  Wannamaker,  1979).  The  jet  enters  the  basin  and  flows 

northeast  to  approximately  4°W,  curves  southward  and  then  splits  (Figure 
1).  Part  of  the  jet  flows  to  the  west  and  is  incorporated  in  an 

antlcyclonlc  gyre,  while  the  remainder  flows  southeast  to  Cape  Tres 
Forcas  and  then  along  the  African  coast  forming  the  southern  periphery 
of  a  cyclonic  circulation. 

Satellite  infrared  imagery  indicates  variations  in  the  shape, 
location  and  intensity  of  the  persistent  anti  cyclonic  gyre  which 

dominates  the  western  Alboran  Sea.  Hydrographic  data  and  satellite 
infrared  images  support  the  year-round  persistence  of  the  gyre,  although 
its  size  and  location  varies  (Stevenson,  1977;  Cheney,  1978;  Wannamaker, 
1979;  Burkov  et  al.,  1979;  Gallagher  et  al . ,  1981).  In  the  eastern 
portion  of  the  Alboran  Sea  a  general  pattern  of  alternating  cyclonic  and 
antlcyclonlc  circulations  has  been  observed  (Cheney,  1978;  Lanoix, 
1974).  Satellite  Imagery  shows  that,  compared  to  the  western  Alboran, 
this  circulation  pattern  is  far  more  variable  and  of  smaller  scale. 

,  A  study  of  circulation  in  the  Alboran  Sea  begins  by  using  the 
simplest  model  capable  of  simulating  major  features  of  the 
circulation.  This  is  a  reduced  gravity  model  in  a  semi -enclosed, 
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rectangular  domain.  It  is  essentially  a  model  of  the  first  baroclinic 
mode  and  does  not  permit  baroclinic  instability  or  the  inclusion  of 
bottom  topography.  We  then  progress  to  the  two  layer  model  eventually 
including  topography  and  the  geometry  of  the  Alboran  Sea.  The  model 

w  * 

formulation  and  parameters  are  discussed, in  Section  2.  In  Sectioa  3  the 
model  results  are  presented  in  terms  of  a  pivotal  experiment  and  some 
deviations  from  it.  Section  3.1  discusses  the  pivotal  experiment.  ''The 
subsections  which  follow  investigate  the  influences  of  (3T/2)  inflow 
angle,  (3.3)  Dort  location,  (3.4)  boundary  effects,  (3.5)  inflow 
vorticity  and  (3.6)  the  addition  of  lower  layer  flow. 


2. 


THE  MOOEL 

A  nonlinear  reduced  gravity  model,  developed  for  the  Gulf  of 
Mexico  by  Hurlburt  and  Thompson  (1980),  has  been  adapted  for  the  Alboran 
Sea.  The  model  equations  are  solved  numerically  using  an  economical 
semi -implicit  scheme.  The  model  consists  of  an  active  upper  layer  and  a 
lower  layer  which  is  infinitely  deep  and  at  rest.  It  is  stably 

stratified  and  has  a  fixed  density  contrast  between  two  immiscible 

layers.  The  model  assumes  a  hydrostatic,  Boussinesq  fluid  in  a  rotating 
right-handed  coordinate  system  on  a  8 -plane.  The  vertically  integrated 
model  equations  are 

3V  ^ 

+  (7  •  Vj  +  Vj  •  7)  Vj  +  k  x  fVj  =  -  g'h1'7h1  + 

(T1  -  ?  )/p  +  AV2V1 

3h1  „  - 
— +  7  •  V1  - 0 

where 

7-|-  i  +  f-  j 

3x  dy  J 

-*•  -*■  *  * 
vi  =  hivi  =  hi  ^ui '  +  vii) 

g'  =  g(P2  -  o})/p 
f  =  fo  +  3(y  -  yQ) 


and  x  and  y  are  tangent-plane  Cartesian  coordinates  with  x  directed 
eastward  and  y  northward,  u^  and  v^  are  the  eastward  and  northward 
velocity  components  in  the  upper  layer,  h1  is  the  upper  layer  thickness, 
t  is  time,  g  is  the  acceleration  due  to  gravity,  ^  is  the  density  of 


seawater  in  layer  i,  fQ  and  y0  are  the  values  of  the  Coriolis  parameter 
and  the  y-coordiante  at  the  southern  boundary,  t 1  is  the  wind  stress, 
and  i2  is  the  interfacial  stress.  The  remaining  parameters  are  defined 
in  Table  1. 


Table  1 

MODEL  PARAMETERS  FOR  THE  PIVOTAL  EXPERIMENT 


Parameter 

Definition 

Value 

A 

eddy  viscosity 

250  m2  sec'1 

8 

(df/dy) 

2  x  lO'^m^sec"1 

fo 

Coriolis  parameter 

8  x  10"5sec'1 

g‘ 

reduced  gravity  due  to  strati¬ 
fication 

.02  m  sec'2 

Hi 

undisturbed  upper  layer  depth 

200  m 

h2 

undisturbed  lower  layer  depth 

00 

LxX 

Ly 

horizontal  dimensions  of  the 
model  domain 

600  km  x  160  km 

X 

X 

Ay 

horizontal  grid  spacing  for  each 
dependent  variable 

10  km  x  5  km 

At 

time  step 

1  hour 

vlin 

inflow  velocity 

100  cm/sec 

a 

angle  of  inflow 

21°  north  of  east 

*s 

inflow  spin  up  time  constant 

30  days 

Figure  1  shows  the  model  domain  superimposed  on  a  map  of  the 
Alboran  Sea.  This  domain  is  600  km  x  160  km  with  10  km  x  5  km  grid 
resolution  for  each  dependent  variable.  Forcing  is  due  solely  to 
prescribed  inflow  through  the  western  port  (Strait  of  Gibraltar). 
Inflow  is  exactly  compensated  by  outflow  through  an  open  eastern 
boundary.  This  is  accomplished  by  allowing  normal  flow  at  the  eastern 
boundary  to  be  self-determined  and  by  imposing  an  integral  constraint  on 
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total  mass  outflow  (Hurlburt  and  Thompson,  1980).  Except  at  the  inflow 
and  outflow  ports,  the  boundaries  are  rigid  and  a  no-slip  condition  is 
prescribed  on  the  tangential  flow.  Along  the  eastern  boundary  the 
tangential  component  is  set  to  zero  one-half  grid  distance  outside  the 
physical  domain. 

The  model  parameters  for  the  pivotal  experiment  are  given  in 
Table  1.  In  this  experiment  the  western  (inflow)  port  is  centered  102 
km  from  the  southern  boundary  and  is  15  km  wide,  a  width  appropriate  for 
the  Strait  of  Gibraltar  at  a  depth  of  100  m.  The  specification  of  the 
model  forcing  is  accomplished  in  either  of  two  ways: 

•  by  prescribing  velocity  (v^)  or 

•  by  prescribing  transport  and  allowing  the  model  to 

determine  the  inflow  velocities. 

The  former  is  used  in  the  pivotal  experiment.  The  total  inflow 

c  o 

transport  used  in  the  model  (2.5  x  10  m  /sec)  is  higher  than 
observed.  This  value  is  necessary  to  drive  a  uni  from  inflow  profile  for 
vx  or  Yj  with  speeds  of  -  100  cm/sec,  given  the  port  width  and  upper 
layer  depth  in  Table  1.  The  inflow  is  spun  up  with  a  time  constant  of 
30  days  to  minimize  the  excitation  of  high  frequency  waves.  The  angle 
of  inflow  was  varied  based  on  direct  observations  (Lacombe,  1961)  and  on 
inferences  from  satellite  imagery.  The  standard  inflow  angle  was  chosen 
to  be  21°  north  of  east  based  on  the  geometric  orientation  of  the  Strait 
of  Gibraltar.  Possibly  important  wind  forcing  (Ovchinnikow  et  al., 
1976;  Mommsen,  1978)  is  neglected  to  allow  focus  on  the  circulation 
driven  by  flow  through  the  Strait  of  Gibraltar. 

Substantial  effort  was  made  to  assure  that  unphysical  aspects 
of  the  model  such  as  the  grid  resolution  and  the  time  step  did  not 
significantly  influence  the  model  solution.  The  open  eastern  boundary 
condition  was  a  special  concern  and  one  important  test  of  its  influence 
is  discussed  in  Section  3.4.  The  integral  constraint  on  the  total  mass 
flux  through  the  eastern  boundary  resulted  in  plane  wave  reflection  of 
sufficient  amplitude  to  pose  a  significant  problem.  The  eddy  viscosity 
(A)  chosen  for  the  model  is  the  minimum  value  which  prevents  any  visible 
oscillation  in  the  solution  due  to  this  reflection.  An  eddy  viscosity 
2t  times  smaller  yielded  nearly  the  same  solution  except  for  some 
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unphysical  oscillations.  A  viscous  boundary  layer  using  a  linear 
interfacial  stress  was  also  applied  near  the  open  eastern  boundary  in  an 
effort  to  damp  the  oscillations  due  to  the  integral  constraint.  The 

O  1 

maximum  value  for  the  drag  coefficient  was  10  sec  at  the  eastern 
boundary.  It  decreased  exponentially  away  from  the  boundary  with  an  e- 
folding  width  of  50  km.  This  aided  only  slightly  in  damping  the 
reflection  from  the  integral  constraint.  Except  for  the  viscous 
boundary  layer,  the  interfacial  stress  was  zero. 

• 


» 


i 


i 


i 
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3.  MODEL  RESULTS 

Over  forty  numerical  experiments  were  performed  for  the  Alboran 
Sea.  Some  preliminary  results  will  be  presented  in  terms  of  a  pivotal 
experiment  and  selected  variations  from  it.  Most  of  the  numerical 
solutions  evolved  to  a  steady  state  in  about  one  year. 

3.1  THE  PIVOTAL  EXPERIMENT 

The  pivotal  experiment  uses  the  parameters  in  Table  1.  Figure 
2  shows  the  steady  state  model  solution  (day  500)  in  terms  of  the 
pycnocline  anomaly  (PA).  The  PA  is  the  divation  of  the  interface 
between  the  upper  and  lower  layers  from  its  initial  flat  position  at  2on 
m  depth.  Downward  deviations  are  positive  (upper  layer  thicker  ti 
initially).  The  most  striking  features  are 

•  a  meandering  current  which  traverses  the  model  domain  fr 
west  to  east, 

•  a  strong  anticyclonic  gyre  in  the  western  249  km  of  the 
basin,  and 

•  a  weak  cyclonic  circulation  to  the  east. 

This  pattern  is  very  similar  to  the  temperature  field  shown  in  a 
satellite  image  and  to  Lanoix's  dynamic  topography  (Figure  1). 


500.  Inflow  angle  is  20°  north  of  east.  Solid  contours  are  positive 
(downward)  deviations.  Dashed  contours  are  negative  (upward) 
deviations.  Contour  interval  is  10  m. 

In  rotating  tank  experiment,  Whitehead  and  Miller  (1979) 
attempted  to  simulate  the  Alboran  Gyre  using  a  density  driven  current. 
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They  suggest  that  the  dimensions  of  the  gyre  depend  on  a  coastline 
feature.  Cape  Tres  Forcas.  It  has  also  been  suggested  by  Porter  (1976) 
that  the  dimensions  of  the  gyre  are  directly  related  to  the  bottom 
topography  with  Alboran  Island  acting  as  the  eastern  boundary  of  the 
gyre.  Yet  the  reduced  gravity  model  is  able  to  simulate  an  Alboran  Gyre 
with  realistic  dimensions,  while  including  neither  coastline 
irregularities  nor  bottom  topography.  The  model  gyre  is  also  a 
persistent  rather  than  a  transient  feature  of  the  flow,  again  in  accord 
with  observations  noted  earlier. 

3.2  THE  EFFECT  OF  INFLOW  ANGLE 

The  circulation  pattern  in  the  pivotal  experiment  (Figure  2) 
and  satellite  imagery  suggest  that  the  inflow  angle  may  exert  an 
important  influence  on  the  meandering  current  and  the  Alboran  Gyre. 
Thus,  a  number  of  experiments  were  carried  out  varying  the  inflow 
angle.  One  experiment  used  the  standard  parameters  from  Table  1  except 
that  the  inflow  entered  the  model  domain  normal  to  the  western 
boundary.  In  the  steady  state  solution  for  this  case  (Figure  3),  the 
jet  enters  the  basin  flowing  due  east,  but  quickly  veers  southward.  The 
anticyclonic  gyre  in  the  western  part  of  the  basin  is  restricted  to  a 
much  smaller  north-south  extent  than  in  the  standard  case.  The  cyclonic 
circulation  to  the  east  is  intensified  and  a  new  anticyclonic  flow 
appears  in  the  eastern  200  km.  This  circulation  pattern  is  similar  to 
that  of  the  sea  surface  temperature  seen  in  the  satellite  infrared 
imagery. 


Figure  3.  PA  of  the  pivotal  case  at  day  500  but  with  due  east  inflow. 
Contour  interval  is  10  m. 
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3.3  PORT  LOCATION 

The  next  set  of  experiments  was  designed  to  observe  the 
importance  of  the  north-south  location  of  the  inflow  port.  If  the 
Strait  of  Gibraltar  was  located  south  of  the  basin  center,  how  would  the 
circulation  pattern  be  affected?  Two  such  experiments  used  the  standard 
parameters  of  Table  1  except  that  the  port  was  centerd  7  km  south  of  the 
center  of  the  western  boudary  and  the  inflow  angles  were  those  of  Figure 
2  (standard  21°  north  of  east)  and  Figure  3  (0°).  The  steady  state 
solutions  are  presented  in  Figure  4.  Figure  4a,  with  the  angles  inflow 
of  the  standard  case,  shows 

•  an  anticyclonic  gyre  smaller  than  that  of  the  standard  case 
in  the  western  part  of  the  basin, 

•  a  strong  cyclonic  circulation  in  the  central  part,  and 

•  a  weak  anticyclonic  circulation  in  the  eastern  part. 

When  the  inflow  enters  flowing  due  east  (Figure  4b)  the  western  gyre  is 
even  smaller  than  in  Figure  4a  and  the  two  gyres  to  the  east  are 
slightly  stronger.  Clearly,  the  entrance  of  the  Atlantic  water  in  the 
northern  half  of  the  Alboran  Sea  facilitates  the  development  of  an 
Alboran  Gyre  of  large  north-south  extent. 

3.4  BOUNDARY  EFFECTS 

The  influences  of  the  domain  size  and  the  open  eastern  boundary 
were  also  investigated.  As  described  in  Section  2,  the  pivotal 
experiment  includes  a  viscous  boundary  layer  near  the  open  eastern 
boundary.  One  test  compared  the  pivotal  experiment  with  this  boundary 
layer  (Figure  2)  to  an  identical  experiment  without  it  (not  shown). 
Except  in  the  region  of  the  viscous  boundary  layer,  the  results  were 
almost  identical.  In  the  experiment  without  the  boundary  layer,  the 
current  maintained  its  cross-sectional  sturcture  as  it  approached  and 
passed  through  the  open  eastern  boundary.  When  it  was  included,  the  PA 
contours  spread  out  in  the  viscous  boundary  layer  and  the  jet  structure 
disintegrated  to  a  more  uniform  flow  (see  Figure  2). 


(KM) 


(b)  (KM)  fey. 


(KM) 


Figure  4.  PA  of  the  pivotal  case  at  day  500  but  with  the  port  centered 
7  km  south  of  the  center  of  the  western  boundary,  (a)  Inflow  angles  21° 
north  of  east;  (b)  due  east  inflow.  Contour  interval  is  10  m. 


Numerical  experiments  were  also  performed  to  determine  if  the 
open  eastern  boundary  was  seriously  distorting  the  solution.  Figure  5 
shows  the  results  of  a  critical  test.  It  compares  two  solutions  which 
differ  only  in  the  east-west  extent  of  the  model  domain.  In  each  case 
there  is  a  viscous  boundary  layer  near  the  open  eastern  boundary.  In 
this  test,  changing  the  location  of  the  open  eastern  boundary  caused 
only  minor  changes  in  the  solution  in  the  western  400  km  of  the  model 
domain.  The  eastern  100  km  in  Figure  5b  differs  from  the  same  region  in 
Figure  5a,  due  mostly  to  the  viscous  boundary  layer  in  the  vicinity  of 
the  open  boundary. 


Figure  5.  PA  of  the  pivotal  case  at  day  190  *ut  with  different 
dimensions,  (a)  180  km  x  800  km  basin;  (b)  180  km  x  500  km  basin.  Scale 
factor  for  these  figures  x:y  is  1:2.5.  Contour  Interval  is  10  m. 


The  effect  of  the  north-south  extent  of  the  basin  was  examined 
by  doubling  the  y-dimension  of  the  standard  experiment  (Figure  2)  and 
keeping  the  port  location  slightly  north  of  the  basin  center.  Figure  6 
shows  the  flow  entering  the  basin  at  the  standard  angel  of  21°  and  then 
curving  southward  in  a  manner  similar  to  Figure  4a.  Even  though  the 
north-south  extent  of  the  basin  has  been  doubled,  the  northern  and 
southern  boundaries  of  the  domain  still  limit  the  north-south  scale  of 
the  current  meanders.  Despite  a  large  increase  in  the  amplitude  of  the 
meanders  in  Figure  6,  the  wavelengths  in  Figure  4a  and  Figure  6  are 
almost  the  same.  A  striking  feature  in  Figure  6  is  the  downstream 
amplification  of  the  current  meanders.  Less  dramatic  examples  of  this 
appear  in  some  of  the  other  figures.  It  should  be  noted  that  this  and 
all  the  other  solutions  are  steady  and  not  unstable. 


Figure  6.  PA  at  day  600  of  an  experiment  where  the  north-south 
dimension  of  the  basin  has  been  extended  to  320  km,  contour  interval  Is 
10  m. 


3.5  SHEAR  AT  INFLOW 

In  all  the  experiments  discussed  so  far,  the  inflow  has  been 
prescribed  as  a  velocity  profile  (v^).  However,  in  this  model  the 
Inflow  may  also  be  prescribed  in  terms  of  transport  ( ) .  In  the  latter 
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case,  the  model  partially  controls  the  inflow  velocity  profile  through 
the  geostrophic  tilt  In  the  interface.  Figure  7  shows  a  steady  state 
solution  for  an  experiment  similar  to  that  shown  in  Figure  4b,  except 
that  is  prescribed  instead  of  Vj.  In  both  cases  the  inflow  is 

eastward.  In  Figure  7  the  prescribed  inflow  transport  is  2.5  x  10^ 
m-Vsec.  This  yields  inflow  velocities  similar  to  Figure  4b,  but  the 
geostrophic  tilt  in  the  interface  introduces  a  shear  (3u/3y)  of  1.54  x 
10"5  sec"*  at  inflow.  Without  the  shear  (Figure  4b)  the  current  turns 
southward  after  inflow,  but  with  the  shear  (Figure  7)  the  current  turns 
northward.  In  addition,  Figure  7  shows  a  larger  Alboran  Gyre  which  is 
further  to  the  east  and  a  small  cyclonic  gyre  in  the  northwest 
corner.  The  possibility  that  vorticity  at  inflow  turns  the  incoming 
Atlantic  water  northward  has  been  suggested  by  Nof  (1978). 


I  KM) 


...  *  v'V" 
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Figure  7.  PA  at  day  500  of  a  case  identical  to  the  case  represented  in 
Figure  4b  except  that  the  model  is  forced  with  a  prescribed  transport. 
Contour  interval  is  10  m. 


3.6  ADDITION  OF  LOWER  LAYER  FLOW 

In  the  Alboran  Sea  region  of  the  Mediterranean,  circulation  may 
be  divided  into  two  layers;  an  upper  eastward  moving  layer,  200  m  deep 
and  a  lower  westward  moving  layer  2400  m  deep.  Figure  8  is  a 
representation  of  the  two  layer  model.  The  equations  used  are  those 
presented  in  Section  2,  where  the  subscript  1  is  replaced  by  2  and  the 
pressure  term  is  replaced  by 
-  h,  v  p, 
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Figure  9a.  Upper  Layer 


Figure  9b.  Lower  Layer 


In  this  first  two  layer  experiment  the  model  dimensions  and 
configuration  are  identical  to  that  of  the  reduced  gravity  model  except 
that  a  uniform  2400  m  deep  lower  layer  has  been  added.  Forcing  in  the 
upper  layer  remains  the  same  as  in  the  reduced  gravity  pivotal  case  ex¬ 
cept  for  the  port  location  which  has  been  moved  to  the  exact  center  of 
the  basin.  The  deep  lower  layer  is  forced  with  a  westward  moving  flow 
from  the  open  eastern  boundary.  This  flow  is  kept  at  a  uniform  .2  cm/ 
sec  across  the  boundary,  a  value  which  agrees  with  observations  that 
near  the  Straits,  lower  layer  flow  Is  as  high  as  10  cm/sec.  Outflow  is 
allowed  only  through  the  narrow  western  port  (Figure  9). 

Figures  10a  and  10b  show  a  comparison  after  300  days  of  model 
integration  of  the  reduced  gravity  (Figure  10a)  and  two  layer  (Figure 
10b)  cases.  Note  there  Is  very  little  difference  in  the  location,  size 
and  shape  of  the  gyre  in  the  two  cases.  Thus,  the  simple  one  layer 
model  seems  to  represent  all  of  the  Important  features  of  the  two  layer 
flat  bottom  case. 


Figure  10a.  Reduced  gravity  model  solution. 


4.  SUMMARY  AND  FUTURE  WORK 

The  nonlinear,  semi -implicit,  reduced  gravity  and  two  layer 
numerical  model  (Hurl hurt  and  Thompson,  1980)  have  been  adapted  to  study 
the  circulation  in  the  Alboran  Sea.  Model  results  using  an  idealized 
rectangular  geometry  (600  km  x  160  km),  no  topography  and  a 
northeastward  inflow  through  the  Strait  of  Gibraltar  show  an 
anitcyclonic  gyre  similar  in  size,  shape  and  location  to  the  Alboran 
Gyre  (Figure  2).  These  results  closely  resemble  the  dynamic  height 
contours  of  Lanoix  (Figure  1)  and  suggest  that  topography  and  particular 
coastline  features  are  not  necessary  to  create  a  gyre  with  realistic 
dimensions  and  location.  However,  model  results  for  the  eastern  Alboran 
show  a  series  of  cyclonic  and  anticyclonic  circulations  of  much  larger 
scale  and  smaller  variability  than  observed.  Additional  numerical 
experiments  showed  the  importance  of  the  inflow  angle  and  inflow 
vorticity  in  determining  the  size  and  location  of  the  Alboran  Gyre. 
When  an  additional  uniformly  flat  lower  layer  with  weak  westward  flow  is 
included,  little  difference  is  observed  between  the  reduced  gravity  and 
two  layer  solutions. 

This  paper  has  presented  preliminary  results  of  an  attempt  to 
model  the  Alboran  Sea.  Future  work  will  include  an  investigation  of  the 
model  dynamics,  more  realistic  models,  in  particular  two  layer  models 
including  geometry  and  topography,  and  interaction  with  a  field 
experiment. 
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